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I SUMMARY 


This paper presents the derivation of the differential stiffness for 
triangular solid of revolution elements. The derivation takes into account the 
! element rigid body rotation only, the rotation being about the circumferential 
t axis. Internal pressurization of a pneumatic tire is used to illustrate the 
(application of this feature. 


INTRODUCTION 


i In the NASTRAN computer code, the rigid format No. 4 is the so-called static 
analysis with differential stiffness. This method of analysis is a first approx- 
imation to large deflection effects. The differential stiffness refers to the 
geometry stiffness which is generated using the stresses in the element which 
are computed from the initial linear stress analysis. The differential stiff- 
ness is available for rods, beams, shear panels, plates, and conical shell ele- 
ments. However, sometimes one desires that the differential stiffness be avail- 
able to solid elements. This is particularly true when one deals with a 
structure having thick wall yet flexible in nature. With this in mind, the 

differential stiffness for a triangular solid-of-revolution element is developed 
and incorporated into NASTRAN. 

In the case of a ring element loaded axi symmetrically, there exists only 
lone rotation of the element which is the rotation about the circumferential axis. 
We will assume that the stress and circumfeeential rotation are uniform through 
the element. Furthermore, in computing the work done by the forces in at n 
element we will neglect the work done in conjunction with the normal and shear 
strains and the work done by the circumferential stress and strain. This is in 
accord with the derivation of the differential stiffness for a triangular 
membrane element presented in reference 1. 
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DERIVATION OF DIFFERENTIAL STIFFNESS MATRIX 


The triangular ring element in NASTRAN assumes the following displacement j 
field (Fig. 1) 

u = $1 + $2' + $ 3 Z (i* a ) 

w = 64 + S 5 r + e 6 z (1-b) 1 


where u, w are the displacements in the direction of the r, 
respectively. Inverting the above yields 

- tv {q} 

where {q} T = w 1 u 2 w 2 u 3 w 3 J 

{3} = L e 1 e 2 e 3 e 4 e 5 e 6 -I 


z coordinates 


(2) 


Due to symmetry the only rotation of the element is the circumferential rotationl 
u» 0 which according to reference 1 is given as 
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Substituting 3_ and from equation (2) into equation (3) we obtain 
5 o 
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and 



r . 
J 


z . 
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4 = r l2 Z 13 ' r l 3 Z 12 

Considering only the circumferential rigid body rotation and neglecting the 
effects of the normal and shear strains, then the work done by forces in an 
element of volume V reduces to 
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(5) 


W = 


\ “e (<T r + ° z ) 


/here a r and o z are the element normal stresses in the r and z directions, 
'espectively. Hence the matrix of the differential stiffness coefficient for 
:hat element is simply 


i 


[K d ] 

*• 0X0 



0 0 
0 0 


( 6 ) 


We then have the differential stiffness matrix in terms of the generalized coor- 
dinate {$} 


[K d ] = [C] T [K d ][C] 

ee too) 


(7) 


The above differential stiffness matrix is then transformed to grid point coor- 
dinates by means of equation (2) as 


[K d ] = [r ] T [K d ] [r. ] (8) 

eq ee Jl 3q J v ' 

This geometric stiffness matrix is then added to the elastic stiffness matrix to 
!yield the total stiffness matrix of the structure. 


APPLICATION TO PNEUMATIC TIRE INFLATION 


j 

We will apply the above formulation to the problem of inflating a pneumatic 
tire. This is an example which can be considered as an axi symmetric problem. 

Fig. 2 shows one-half of a typical meridian cross-section of a bias pneumatic 
tire. Its structural elements consist of (1) rubber which spreads from the 
itread area to shoulder area and to sidewall area and to the bead area, (2) casing 
or carcass which is the main load carrying member, and (3) bead, which seats 
the tire to the wheel rim. The carcass is a multi-layer composite which consists 
of a number of plies stacked together The ply is a unidirectional composite 
made of many flexible high modulus cord of natural textile, synthetic polymer, 
[glass fiber embedded in and bonded to a matrix of low-modulus polymeric material 
'such as rubber. The plies are usually stacking together in a pair of symmetric 
bias cord angle with respect to the circumferential axis. 

j If we know the properties of the matrix and cord materials and their volume 
ratios, then the properties of the composite ply in its principal directions, 
directions 1 and 2 in Fig. 3, can be determined by the well-known Halpin-Tsai 
equation, reference 2. The ply is assumed to be transversely isotropic, i.e. 
Young's modulus E 33 = E 22 and Poisson's ratios v 13 =v 12 , v =v 23 . The well-known 

transformation matrix for elastic constants through a rotation is then applied 
to obtain the ply stiffness matrix in the laminate coordinate system x, y. For 
the present case, they are the meridional and circumferential directions. The 
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cord angle a which is also the angle of rotation of the transformation is cal- 
culated from the following equation known as "lift equation" in the tire 
industry. | 

r 1 = cos S 

r c 1 + e cos a 

where a = cord angle in the body of tire which takes the molded shape 
g = cord angle in the initial cylinder before molding 
r = radius coordinate of the location in the body of tire 
r = radius of the initial cylinder before molding 

e = estimated cord strain i 

The total effective carcass stiffness matrix is then obtained by the method in ] 
reference 3. In general, the laminate can be satisfactorily represented by two-] 
ply laminate of (+a) and (-a) cord orientations. j 

The material properties used in our example are 

cord modulus E c = 0. 52376x1 0 5 kg/cm 2 (7. 45x1 0 5 psi) 

cord Poisson's ratio v a 0.3 

rubber modulus E r = 52.7278 kg/cm 2 (750 psi) 

rubber Poisson's ratio v f = 0.48 

Other data used are 

g = 57 degrees 

r c = 21.4846 cm (8.4585 in.) 

The cord strain e is assumed to be constant and its value at the crown is used. 
The radius of the ply center line at the crown is 32.8308 cm (12.9255 in.) and 
the final cord angle at that point is assumed to be 32 degrees with reference 
to circumferential circle. The tire is subjected to internal pressure of 

1.6873 kg/cm 2 (24 psi). For simplicity, the tire is assumed to be fixed at a 
section close to the bead. In all, 488 elements and 317 grid points are used 
for one-half of the tire cross-section. Figure 4 shows the shape of a 
deformed and undeformed tire section for the linear solution and for the 
differential stiffness solution. The radial displacement of the outer ply line 
using NASTRAN is 1.3785 cm (0.54277 in.) compared with 1.2217 cm (0.481 in.) 
from a shell theory. The shell theory neglects the rubber in the tread, 
shoulder, sidewall, and bead and only considers the carcass as a thin elastic 
shell of orthotropic material. The increase in the tire outer radius due to 
inflation is computed by NASTRAN to be 1.3419 cm (0.52829 in.) as compared with 
1.0976 cm (0.448 in.) from measurement. The radial displacement of the same 
point from NASTRAN linear solution is 1.7499 cm (0.68895 in.). The inflated 
dimension of the outer carcass line in the z direction is 10.3045 cm (4.0569 
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yth f nf m tho S I?J N ^II d 9,9847 c "]. (3-931 in.) from shell theory. The section 
Pth of the tire(the maximum dimension of the tire outer profile in the z 

t h C 20°?7 7 ? r r « n f« 1 0 ? I s 21 V 116 cm (8 ' 3134 in -> from NASTRAN as compared 

the differential sti4fL^°rPQnit me ^ SUrement i - As can be seen> the i nc1u sion 
tne differential stiffness results in a great improvement over the linear 

lut,on. However there are still discrepancies in the various results The 

screpancies can be attributed to factors such as material characterization 

sumptions on the boundary conditions, initial shape, cord anqles cord strain 

d the possibility of the residual stresses in the a^Sal tire * 


CONCLUDING REMARKS 


r trLnau? a rcoi?5 e of here ? h J- derivation of the differential stiffness matrix 

inJ the iinJi ™ £ l mJM 1 - 0 " elen,ents * J ts application is illustrated 
ing the inflation of a bias tire as an example. The differential stiffness 

a ring eTement is particularly useful when dealing with structure which has 

1 evamli/nf ri ex ] ble *: n Si 9 !?. to undergo large displacement and rotation. 
F®J am P le ° f the tore inflation demonstrates this point. The numerical 
kults in the present example are quite satisfactory. 


As mentioned earlier, the derivation of the differential stiffness present- 
iin this paper is based on the information contained in the manual for Level 

1 reference 1, since it was available to this author at that time. This 
thor has learned since then that Level 15.5 (reference 4) has a refined 
nvation of the differential stiffness which includes the effects of the 

strains in computing the energy. Hence, if a new derivation of 

2 differential stiffness for a triangular ring element were made to include 
2 effects of the normal and shear strains and also to include the effect of 

2 circumferential strain, it would probably give a better solution than the 
esent one. 
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Figure 1.- Triangular solid of revolution element. 



Figure 2.- Meridian cross section of a tire. 
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LINEAR SOLUTION 
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(a) Linear solution. 

Figure 4.- Deformed and undeformed tire 


cross section. 




1 / 

^-DIFFERENTIAL STIFFNESS 
SOLUTION 


(t>) Differential stiffness. 
Figure 4.- Concluded. 
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